Machine learning identifies multi-parametric functional PET/MR imaging cluster to predict radiation resistance in preclinical head and neck cancer models.

PURPOSE: Tumor hypoxia and other microenvironmental factors are key determinants of treatment resistance. Hypoxia positron emission tomography (PET) and functional magnetic resonance imaging (MRI) are established prognostic imaging modalities to identify radiation resistance in head-and-neck cancer (HNC). The aim of this preclinical study was to develop a multi-parametric imaging parameter specifically for focal radiotherapy (RT) dose escalation using HNC xenografts of different radiation sensitivities. METHODS: A total of eight human HNC xenograft models were implanted into 68 immunodeficient mice. Combined PET/MRI using dynamic [18F]-fluoromisonidazole (FMISO) hypoxia PET, diffusion-weighted (DW), and dynamic contrast-enhanced MRI was carried out before and after fractionated RT (10 × 2 Gy). Imaging data were analyzed on voxel-basis using principal component (PC) analysis for dynamic data and apparent diffusion coefficients (ADCs) for DW-MRI. A data- and hypothesis-driven machine learning model was trained to identify clusters of high-risk subvolumes (HRSs) from multi-dimensional (1-5D) pre-clinical imaging data before and after RT. The stratification potential of each 1D to 5D model with respect to radiation sensitivity was evaluated using Cohen's d-score and compared to classical features such as mean/peak/maximum standardized uptake values (SUVmean/peak/max) and tumor-to-muscle-ratios (TMRpeak/max) as well as minimum/valley/maximum/mean ADC. RESULTS: Complete 5D imaging data were available for 42 animals. The final preclinical model for HRS identification at baseline yielding the highest stratification potential was defined in 3D imaging space based on ADC and two FMISO PCs ([Formula: see text]). In 1D imaging space, only clusters of ADC revealed significant stratification potential ([Formula: see text]). Among all classical features, only ADCvalley showed significant correlation to radiation resistance ([Formula: see text]). After 2 weeks of RT, FMISO_c1 showed significant correlation to radiation resistance ([Formula: see text]). CONCLUSION: A quantitative imaging metric was described in a preclinical study indicating that radiation-resistant subvolumes in HNC may be detected by clusters of ADC and FMISO using combined PET/MRI which are potential targets for future functional image-guided RT dose-painting approaches and require clinical validation.

Comparison of patient stratification by computed tomography radiomics and hypoxia positron emission tomography in head-and-neck cancer radiotherapy.

BACKGROUND AND PURPOSE: Hypoxia Positron-Emission-Tomography (PET) as well as Computed Tomography (CT) radiomics have been shown to be prognostic for radiotherapy outcome. Here, we investigate the stratification potential of CT-radiomics in head and neck cancer (HNC) patients and test if CT-radiomics is a surrogate predictor for hypoxia as identified by PET. MATERIALS AND METHODS: Two independent cohorts of HNC patients were used for model development and validation, HN1 (n = 149) and HN2 (n = 47). The training set HN1 consisted of native planning CT data whereas for the validation cohort HN2 also hypoxia PET/CT data was acquired using [18F]-Fluoromisonidazole (FMISO). Machine learning algorithms including feature engineering and classifier selection were trained for two-year loco-regional control (LRC) to create optimal CT-radiomics signatures.Secondly, a pre-defined [18F]FMISO-PET tumour-to-muscle-ratio (TMRpeak ≥ 1.6) was used for LRC prediction. Comparison between risk groups identified by CT-radiomics or [18F]FMISO-PET was performed using area-under-the-curve (AUC) and Kaplan-Meier analysis including log-rank test. RESULTS: The best performing CT-radiomics signature included two features with nearest-neighbour classification (AUC = 0.76 ± 0.09), whereas AUC was 0.59 for external validation. In contrast, [18F]FMISO TMRpeak reached an AUC of 0.66 in HN2. Kaplan-Meier analysis of the independent validation cohort HN2 did not confirm the prognostic value of CT-radiomics (p = 0.18), whereas for [18F]FMISO-PET significant differences were observed (p = 0.02). CONCLUSIONS: No direct correlation of patient stratification using [18F]FMISO-PET or CT-radiomics was found in this study. Risk groups identified by CT-radiomics or hypoxia PET showed only poor overlap. Direct assessment of tumour hypoxia using PET seems to be more powerful to stratify HNC patients.

PET/MRI and genetic intrapatient heterogeneity in head and neck cancers.

PURPOSE: The relation between functional imaging and intrapatient genetic heterogeneity remains poorly understood. The aim of our study was to investigate spatial sampling and functional imaging by FDG-PET/MRI to describe intrapatient tumour heterogeneity. METHODS: Six patients with oropharyngeal cancer were included in this pilot study. Two tumour samples per patient were taken and sequenced by next-generation sequencing covering 327 genes relevant in head and neck cancer. Corresponding regions were delineated on pretherapeutic FDG-PET/MRI images to extract apparent diffusion coefficients and standardized uptake values. RESULTS: Samples were collected within the primary tumour (n = 3), within the primary tumour and the involved lymph node (n = 2) as well as within two independent primary tumours (n = 1). Genetic heterogeneity of the primary tumours was limited and most driver gene mutations were found ubiquitously. Slightly increasing heterogeneity was found between primary tumours and lymph node metastases. One private predicted driver mutation within a primary tumour and one in a lymph node were found. However, the two independent primary tumours did not show any shared mutations in spite of a clinically suspected field cancerosis. No conclusive correlation between genetic heterogeneity and heterogeneity of PET/MRI-derived parameters was observed. CONCLUSION: Our limited data suggest that single sampling might be sufficient in some patients with oropharyngeal cancer. However, few driver mutations might be missed and, if feasible, spatial sampling should be considered. In two independent primary tumours, both lesions should be sequenced. Our data with a limited number of patients do not support the concept that multiparametric PET/MRI features are useful to guide biopsies for genetic tumour characterization.

Automatic VMAT planning for post-operative prostate cancer cases using particle swarm optimization: A proof of concept study.

OBJECTIVE: To investigate the potential of Particle Swarm Optimization (PSO) for fully automatic VMAT radiotherapy (RT) treatment planning. MATERIAL AND METHODS: In PSO a solution space of planning constraints is searched for the best possible RT plan in an iterative, statistical method, optimizing a population of candidate solutions. To identify the best candidate solution and for final evaluation a plan quality score (PQS), based on dose volume histogram (DVH) parameters, was introduced. Automatic PSO-based RT planning was used for N = 10 postoperative prostate cancer cases, retrospectively taken from our clinical database, with a prescribed dose of EUD = 66 Gy in addition to two constraints for rectum and one for bladder. Resulting PSO-based plans were compared dosimetrically to manually generated VMAT plans. RESULTS: PSO successfully proposed treatment plans comparable to manually optimized ones in 9/10 cases. The median (range) PTV EUD was 65.4 Gy (64.7-66.0) for manual and 65.3 Gy (62.5-65.5) for PSO plans, respectively. However PSO plans achieved significantly lower doses in rectum D2% 67.0 Gy (66.5-67.5) vs. 66.1 Gy (64.7-66.5, p = 0.016). All other evaluated parameters (PTV D98% and D2%, rectum V40Gy and V60Gy, bladder D2% and V60Gy) were comparable in both plans. Manual plans had lower PQS compared to PSO plans with -0.82 (-16.43-1.08) vs. 0.91 (-5.98-6.25). CONCLUSION: PSO allows for fully automatic generation of VMAT plans with plan quality comparable to manually optimized plans. However, before clinical implementation further research is needed concerning further adaptation of PSO-specific parameters and the refinement of the PQS.

Potentials and challenges of diffusion-weighted magnetic resonance imaging in radiotherapy.

PURPOSE: To review the potential and challenges of integrating diffusion weighted magnetic resonance imaging (DWI) into radiotherapy (RT). CONTENT: Details related to image acquisition of DWI for RT purposes are discussed, along with the challenges with respect to geometric accuracy and the robustness of quantitative parameter extraction. An overview of diffusion- and perfusion-related parameters derived from mono- and bi-exponential models is provided, and their role as potential RT biomarkers is discussed. Recent studies demonstrating potential of DWI in different tumor sites such as the head and neck, rectum, cervix, prostate, and brain, are reviewed in detail. CONCLUSION: DWI has shown promise for RT outcome prediction, response assessment, as well as for tumor delineation and characterization in several cancer types. Geometric and quantification robustness is challenging and has to be addressed adequately. Evaluation in larger clinical trials with well designed imaging protocol and advanced analysis models is needed to develop the optimal strategy for integrating DWI in RT.

Assessment of image quality of a radiotherapy-specific hardware solution for PET/MRI in head and neck cancer patients.

BACKGROUND AND PURPOSE: Functional PET/MRI has great potential to improve radiotherapy planning (RTP). However, data integration requires imaging with radiotherapy-specific patient positioning. Here, we investigated the feasibility and image quality of radiotherapy-customized PET/MRI in head-and-neck cancer (HNC) patients using a dedicated hardware setup. MATERIAL AND METHODS: Ten HNC patients were examined with simultaneous PET/MRI before treatment, with radiotherapy and diagnostic scan setup, respectively. We tested feasibility of radiotherapy-specific patient positioning and compared the image quality between both setups by pairwise image analysis of 18F-FDG-PET, T1/T2-weighted and diffusion-weighted MRI. For image quality assessment, similarity measures including average symmetric surface distance (ASSD) of PET and MR-based tumor contours, MR signal-to-noise ratio (SNR) and mean apparent diffusion coefficient (ADC) value were used. RESULTS: PET/MRI in radiotherapy position was feasible - all patients were successfully examined. ASSD (median/range) of PET and MR contours was 0.6 (0.4-1.2) and 0.9 (0.5-1.3) mm, respectively. For T2-weighted MRI, a reduced SNR of -26.2% (-39.0--11.7) was observed with radiotherapy setup. No significant difference in mean ADC was found. CONCLUSIONS: Simultaneous PET/MRI in HNC patients using radiotherapy positioning aids is clinically feasible. Though SNR was reduced, the image quality obtained with a radiotherapy setup meets RTP requirements and the data can thus be used for personalized RTP.

Overlap of highly FDG-avid and FMISO hypoxic tumor subvolumes in patients with head and neck cancer.

BACKGROUND: PET imaging may be used to personalize radiotherapy (RT) by identifying radioresistant tumor subvolumes for RT dose escalation. Using the tracers [18F]-fluorodeoxyglucose (FDG) and [18F]-fluoromisonidazole (FMISO), different aspects of tumor biology can be visualized. FDG depicts various biological aspects, e.g., proliferation, glycolysis and hypoxia, while FMISO is more hypoxia specific. In this study, we analyzed size and overlap of volumes based on the two markers for head-and-neck cancer patients (HNSCC). MATERIAL AND METHODS: Twenty five HNSCC patients underwent a CT scan, as well as FDG and dynamic FMISO PET/CT prior to definitive radio-chemotherapy in a prospective FMISO dose escalation study. Three PET-based subvolumes of the primary tumor (GTVprim) were segmented: a highly FDG-avid volume VFDG, a hypoxic volume on the static FMISO image acquired four hours post tracer injection (VH) and a retention/perfusion volume (VM) using pharmacokinetic modeling of dynamic FMISO data. Absolute volumes, overlaps and distances to agreement (DTA) were evaluated. RESULTS: Sizes of PET-based volumes and the GTVprim are significantly different (GTVprim>VFDG>VH >VM; p < .05). VH is covered by VFDG or DTAs are small (mean coverage 74.4%, mean DTA 1.4 mm). Coverage of VM is less pronounced. With respect to VFDG and VH, the mean coverage is 48.7% and 43.1% and the mean DTA is 5.3 mm and 6.3 mm, respectively. For two patients, DTAs were larger than 2 cm. CONCLUSIONS: Hypoxic subvolumes from static PET imaging are typically covered by or in close proximity to highly FDG-avid subvolumes. Therefore, dose escalation to FDG positive subvolumes should cover the static hypoxic subvolumes in most patients, with the disadvantage of larger volumes, resulting in a higher risk of dose-limiting toxicity. Coverage of subvolumes from dynamic FMISO PET is less pronounced. Further studies are needed to explore the relevance of mismatches in functional imaging.

Comparison of DCE-MRI kinetic parameters and FMISO-PET uptake parameters in head and neck cancer patients.

PURPOSE: Tumor hypoxia is a major cause of radiation resistance, often present in various solid tumors. Dynamic [18 F]-fluoromisonidazole (FMISO) PET imaging is able to reliably assess tumor hypoxia. Comprehensive characterization of tumor microenvironment through FMISO-PET and dynamic contrast enhanced (DCE) MR multimodality imaging might be a valuable alternative to the dynamic FMISO-PET acquisition. The aim of this work was to explore the correlation between the FMISO-PET and DCE-MRI kinetic parameters. METHODS: This study was done on head and neck cancer patients (N = 6), who were imaged dynamically with FMISO-PET and DCE-MRI on the same day. Images were registered and analyzed for kinetics on a voxel basis. FMISO-PET images were analyzed with the two-tissue compartment three rate-constant model. Additionally, tumor-to-muscle ratio (TMR) maps were evaluated. DCE-MRI was analyzed with the extended Tofts model. Voxel-wise Pearson's coefficients were calculated for each patient to assess pairwise parameter correlations. RESULTS: Median correlations between FMISO uptake parameters and DCE-MRI kinetic parameters varied across the parameter pairs in the range from -0.05 to 0.71. The highest median correlation of r = 0.71 was observed for the pair Vb -vp , while the K1 -Ktrans median correlation was r = 0.45. Median correlation coefficients for the K1 -vp and the Ki -Ktrans pairs were r = 0.42 and r = 0.32, respectively. Correlations between FMISO uptake rate parameter Ki and DCE-MRI kinetic parameters varied substantially across the patients, whereas correlations between the FMISO and DCE-MRI vascular parameters were consistently high. Median TMR-K1 and TMR-Ktrans correlations were r = 0.52 and r = 0.46, respectively, but varied substantially across the patients. CONCLUSIONS: Based on this clinical evidence, we can conclude that the vascular fraction parameters obtained through DCE-MRI kinetic analysis or FMISO kinetic analysis measure the same biological property, while other kinetic parameters are unrelated. These results might be useful in the design of future clinical trials involving FMISO-PET/DCE-MR multimodality imaging for the assessment of tumor microenvironment.

Competing salt effects on phase behavior of protein solutions: tailoring of protein interaction by the binding of multivalent ions and charge screening.

The phase behavior of protein solutions is affected by additives such as crowder molecules or salts. In particular, upon addition of multivalent counterions, a reentrant condensation can occur; i.e., protein solutions are stable for low and high multivalent ion concentrations but aggregating at intermediate salt concentrations. The addition of monovalent ions shifts the phase boundaries to higher multivalent ion concentrations. This effect is found to be reflected in the protein interactions, as accessed via small-angle X-ray scattering. Two simulation schemes (a Monte Carlo sampling of the counterion binding configurations using the detailed protein structure and an analytical coarse-grained binding model) reproduce the shifts of the experimental phase boundaries. The results support a consistent picture of the protein interactions responsible for the phase behavior. The repulsive Coulomb interaction is varied by the binding of multivalent counterions and additionally screened by any increase of the ionic strength. The attractive interaction is induced by the binding of multivalent ions, most likely due to ion bridging between protein molecules. The overall picture of these competing interactions provides interesting insight into possible mechanisms for tailoring interactions in solutions via salt effects.

A strategy for multimodal deformable image registration to integrate PET/MR into radiotherapy treatment planning.

BACKGROUND: Combined positron emission tomography (PET)/magnetic resonance imaging (MRI) is highly promising for biologically individualized radiotherapy (RT). Hence, the purpose of this work was to develop an accurate and robust registration strategy to integrate combined PET/MR data into RT treatment planning. Material and methods. Eight patient datasets consisting of an FDG PET/computed tomography (CT) and a subsequently acquired PET/MR of the head and neck (HN) region were available. Registration strategies were developed based on CT and MR data only, whereas the PET components were fused with the resulting deformation field. Following a rigid registration, deformable registration was performed with a transform parametrized by B-splines. Three different optimization metrics were investigated: global mutual information (GMI), GMI combined with a bending energy penalty (BEP) for regularization (GMI+ BEP) and localized mutual information with BEP (LMI+ BEP). Different quantitative registration quality measures were developed, including volumetric overlap and mean distance measures for structures segmented on CT and MR as well as anatomical landmark distances. Moreover, the local registration quality in the tumor region was assessed by the normalized cross correlation (NCC) of the two PET datasets. RESULTS: LMI+ BEP yielded the most robust and accurate registration results. For GMI, GMI+ BEP and LMI+ BEP, mean landmark distances (standard deviations) were 23.9 mm (15.5 mm), 4.8 mm (4.0 mm) and 3.0 mm (1.0 mm), and mean NCC values (standard deviations) were 0.29 (0.29), 0.84 (0.14) and 0.88 (0.06), respectively. CONCLUSION: Accurate and robust multimodal deformable image registration of CT and MR in the HN region can be performed using a B-spline parametrized transform and LMI+ BEP as optimization metric. With this strategy, biologically individualized RT based on combined PET/MRI in terms of dose painting is possible.